Equation of state for shock-compressed porous molybdenum from first-principles 

mean-field potential calculations 
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The Hugoniot curves for shock-compressed molybdenum with initial porosities of 1.0, 1.26, 1.83, 
and 2.31 are theoretically investigated. The method of calculations combines the first-principles 
treatment for zero- and finite-temperature electronic contribution and the mean-field-potential ap- 
proach for the ion-thermal contribution to the total free energy. Our calculated results reproduce 
the Hugoniot properties of porous molybdenum quite well. At low porosity, in particular, the cal- 
culations show a complete agreement with the experimental measurements over the full range of 
data. For the two large porosity values of 1.83 and 2.31, our results are well in accord with the 
experimental data points up to the particle velocity of 3.5 km/s, and tend to overestimate the shock- 
wave velocity and Hugoniot pressure when further increasing the particle velocity. In addition, the 
temperature along the principal Hugoniot is also extensively investigated for porous molybdenum. 
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Molybdenum (Mo) is a high technology metal with 
wide engineering applications for its thermal and me- 
chanical strength, and also chemical resistances. It has 
been used as flyer and/or back reflector for shock wave ex- 
periments to provide accurate standard of high-pressure 
equation of state (EOS)*^. From this view, the precise 
EOS of Mo is critical for theoretical and practical impli- 
cations. To obtain the EOS of Mo, the shock-compressed 
crystal Mo has been investigated by using different ex- 
periment systems^&I&^i^^HiIi. The EOS shock- 
wave data have now been obtained at pressures ranging 
from 0.1 GPa up to a few TPa. Furthermore, to extend 
the EOS to regions of higher internal energy and temper- 
ature, the porous samples of Mo have also been exper- 
imentally investigated by Bakanova et al^ and Trunin 
et al.—. Compared to the extensive theoretical calcula- 
tions and analysis of the shock-compressed crystal Mo, 
the EOS shock-wave data of the porous Mo remain yet 
to be theoretically exploited and understood, which is a 
main driving force for our present study. 

In this paper we calculate the Hugoniots of porous Mo 
with experimentally relevan t 15 ' 16 porosities m=1.0, 1.26, 
1.83, and 2.31 using the mean-field potential (MFP) ap- 
proach facilitated with the first-principles calculations. 
The MFP approach, which will be briefly described be- 
low, was initiated by Wang and co-workers^ 7 -, and proved 
to provide a numerically convenient way to take into ac- 
count the ion-thermal contribution in the total internal 
energy of the metal based on the first-principles calcula- 
tion of zero-temperature internal energy. A lot of elemen- 
tal metals have been successfully teste d 17 i 17 ' 18 i 19 i 20 i 21 i 22 , 
and a general agreement between the MFP calculations 
and the experimental Hugoniots has been achieved un- 
der low shock temperatures (compared to the Fermi tem- 
perature Tp of the valence electrons) and low Hugoniot 
pressures (typically within 0.1 TPa). Also, the shock- 
compressed porous carbon has been investigated using 
this approach^. Our present results of the Hugoniots 
for the porous Mo with different porosities show good 



agreement with the attainable experimental data from 
several groups. 

Now we start with a brief review of the MFP approach. 
For a system with a given averaged atomic volume V and 
temperature T, the Helmholtz free-energy F(V, T) per 
atom can be written as 

F(V,T) = E C (V) + F ion (V, T) + F e i (V, T) , (1) 

where E c represents the 0-K total energy which is 
obtained from ab initio electronic structure calcula- 
tions, F e i is the free energy due to the thermal excita- 
tion of electrons, and Fi on is the ionic vibrational free 
energy which is evaluated from the partition function 
Zion=exY>(—N Fi on /kBT). Here N is the total number of 
lattice ions. In the mean-field approximation, the classi- 
cal Zi on is given bj*2^ 
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The essential of the MFP approach is that the mean-field 
potential <?(r, V) is simply constructed in terms of E c as 
follows^ 7 - 



g(r,V) = ^[E c (R + r) + E c (R 



r) - 2E C (R)} , (3) 



where r represents the distance that the lattice ion 
deviates from its equilibrium position R. It should 
be mentioned that the well-known Dugdale-MacDonald 
expression^ 5 - for the Griineisen parameter can be derived 



by expanding g(r, V) to order r 
formulated as 



Then, Fi on can be 
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with 



v f (V,T) = ^J exp 



9(r,V) 
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r 2 dr. 
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When the electron-phonon interaction and the mag- 
netic contribution are neglected, the electronic contribu- 
tion to the free energy is F e i =E e i — TS e i , where the bare 
electronic entropy S e i takes the form 2 ^ 

S el (V, T) = -k B J n(e, V) [f In / + (1 - /) ln(l - /)] de, 

(6) 

where n(e, V) is the electronic density of states (DOS) 
and / is the Fermi distribution function. With respect 
to Eq. (6), the energy E e \ due to electron excitations can 
be expressed as 



E e i (V T) = J n(e, V)fede - J " n(e, V)ede, 



(7) 



where is the Fermi energy. Given the Helmholtz free- 
energy F(V, T), the other thermodynamic functions such 
as the entropy S — —(dF/dT)v, the internal energy 
E=F + TS, the pressure P=-(dF/dV) T , and the Gibbs 
free energy G=F + PV, can be readily calculated. 

In all of our calculations we take the molybde- 
num structure to be body-centered-cubic (bcc) struc- 
ture {Im3m). The 0-K total energy was calculated us- 
ing the full-potential linearized augmented plane wave 
(LAPW) method 27 in the generalized gradient approxi- 
mation (GGA) 2 -. We used constant muffin-tin radii of 
2.05ao (do is the Bohr radius). The plane wave cutoff 
K cut is determined from R mt x K cut —10.Q. 4000 k points 
in the full zone are used for reciprocal-space integrations. 
The calculations were performed for atomic lattice pa- 
rameter ranging from 4.8ao to 6.8ao. For the highest 
calculated pressure the atomic volume has exceeded the 
touching sphere limit, an extrapolation for getting the 0- 
K total energy points by Morse function has been done. 

The P-V Hugoniot was obtained from the Rankine- 
Hugoniot relation P(V P - V) /2=E - E p for internal en- 
ergy E, pressure P, and volume V, which are achieved by 
shock from initial conditions V P and E p for the present 
porous Mo. Usually, the density of porous materials of 
different porosities (to < 2 ~ 3) becomes at once the 
same as that of crystal under about 1.0 GPa shock- wave 
compression. In this case, it is a good approximation to 
take the initial internal energy Eq for porous Mo to be 
exactly the same as the internal energy Eq for crystal 
Mo. The V P is V p =mVo, where m is the initial porosity 
and Vo is the ambient volume of nonporous single-crystal 
Mo. Note that during the calculations, we have neglected 
shock melting and the phase dependence of the high- 
temperature equation of state, which is in accord with 
the finding of Mitchell et al£2- that the effects of shock 
melting on the P-V Hugoniots of several reference metals 
were too weak to be observed. 

Figure 1 shows the calculated results of U p -U s diagram 
for m=1.0, 1.26, 1.83, and 2.31. Here U p is particle ve- 
locity and U s is shock-wave velocity. For comparison, 
the experimental data£ & 13 ' 16 are also summarized in Fig. 
1. For m=1.0 (crystal Mo), our calculated U p -U s Hugo- 
niot agrees well with the three groups of experimental 
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FIG. 1: (Color online) Calculated shock- wave velocity versus 
particle velocity of molybdenum for initial porosities m—1.0, 
1.26, 1.83, and 2.31. Different calculated curves represent 
different porosities. For comparison, the experimental data 
for m=1.0 (McQueen et a/.-, LASIj£, and Hixson et alM), 
1.26, 1.83, and 2.31 (Trunin et aZ.— ) are also shown in the 
figure. 



data. Of these three groups of experimental data, the 
data points of McQueen et al& are obtained by using 
explosive system. The data points Marsh et alfi (LASL 
Shock Hugoniot Data) are in the same particle velocity 
range from 0.4 to 2.5 km/s. The data points of Hixson 
et alM- for to=1.0 are obtained by using the two-stage 
light-gas gun facility at Los Alamos National Labora- 
tory (LANL) with U p ranging from 2.2 to 4.4 km/s. The 
least-squares fit to these three groups of data gives the 
relation U s =5. 109+1. 247U P at to=1.0, while our calcu- 
lated results gives ?J S =4. 915+1. 347f7 p . For m=1.26, it 
reveals in Fig. 1 that our calculated U s versus U p rela- 
tion agrees with the data points of Trunin et al^ very 
well. For the two larger porosities m=1.83 and to=2.31, 
our results agree well with the data points of Trunin et 
al*£ up to the particle velocity of 3.5 km/s. At even 
more higher particle velocity (U p >4.0 km/s), the scarce 
attainable experimental data for to=1.83 and to=2.31 
show that the amplitude of the shock-wave velocity is 
underestimated by our calculations, and this deviation 
seems to be growing with increasing the particle velocity. 

Figure 2 shows the calculated Hugoniot pressures as 
functions of the particle velocity U p for different porosi- 
ties. The experimental data£ & 13 i 16 are also plotted for 
comparison. For crystal Mo at equilibrium, i.e., to=1.0, 
one can see that the theory agrees well with the experi- 
ments over the full range of data. For to=1.26, our cal- 
culated results are also in excellent agreement with the 
experiment by Trunin et al^. Note that for this porosity, 
there still lacks the data points at high pressure, and our 
calculation above the pressure of 147 GPa needs to be ex- 
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FIG. 2: (Color online) Calculated Hugoniot pressure versus 
particle velocity of molybdenum for initial porosities m—1.0, 
1.26, 1.83, and 2.31. The experimental data points from sev- 
eral groups are also shown in the same way as occurred in 
Fig. 1. 



perimentally verified in the future. For m=1.83 and 2.31, 
our results agree well with the experimental data points 
up to the Hugoniot pressure of 100 GPa. With further 
increase of the particle velocity, a slight deviation of the 
present calculated Hugoniot pressure from the measure- 
ment occurs, with the latter systematically larger than 
the former at given U p . 

Figure 3 shows our results of pressure dependence of 
Hugoniot volume (scaled by ambient atomic volume Vo of 
the crystal Mo) for different porosities. Again, the attain- 
able experimental data*i* i 13 i 16 are also illustrated in Fig. 

3 for comparison. For to=1.0 and 1.26, the agreement be- 
tween our calculated results and experimental measure- 
ment is obviously good. For the porosities m=1.83 and 
2.31, on the other side, the experimental points become 
scattered and are difficult to fit in a smooth curve. For 
this reason, the difference between our calculated V-P 
relation and the experimental data is somewhat enlarged 
compared to the U p -U s and U p -P relations. The scatter- 
ing of data points is most likely due to the fact that it is 
more difficult for larger porosity to prepare the samples 
with the same initial porosity for several experiments. It 
should be mentioned that the V-P relation is more sen- 
sitive to the difference between theoretical results and 
experimental points than that of U p -U s and U p -P rela- 
tions. 

The calculated Hugoniot temperature is shown in Fig. 

4 as a function of pressure P for several values of porosity. 
For comparison, the previous calculations by McQueen 
et aii* and Hixson et ai.— for the case of m=1.0 are also 
shown in Fig. 4. Clearly, our results for m=1.0 is in good 
agreement with the previous results. One can see from 
Fig. 4 that the larger the initial porosity is, the steeper 
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FIG. 3: (Color online) Calculated Hugoniot pressure versus 
relative atomic volume of molybdenum for initial porosities 
m=1.0, 1.26, 1.83, and 2.31. The experimental data points 
from several groups are also shown in the same way as oc- 
curred in Fig. 1. 
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FIG. 4: (Color online) Calculated Hugoniot temperatures as 
functions of Hugoniot pressure P for molybdenum with initial 
porosities m=1.0, 1.26, 1.83, and 2.31. The filled circle and 
the star points represent theoretical results of McQueen et 
ai.— and Hixson et ai.™, respectively. 



the temperature curve is. This is due to that the porous 
material can absorb and transform to heat the energy of 
shock- wave in the process of shock compression. 

Finally, figure 5 shows the calculated Hugoniot tem- 
perature as a function of atomic volume for the values of 
m=1.00, 1.26, 1.83, and 2.31. The previous theoretical 
results*** for m=1.0 are also plotted in Fig. 5 for com- 
parison. The agreement between our results and those of 
the two group s 4 ' 13 is very good. Due to the reason men- 
tioned above, the Hugoniot temperature increases rapidly 
with the increase of initial porosity at a given relative 
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FIG. 5: (Color online) Hugoniot temperature versus rela- 
tive volume of molybdenum for initial porosities m—1.0, 1.26, 
1.83, and 2.31. The filled circle and star points represent 
theoretical results of McQueen et aZ.— and Hixson et aZ.— , 
respectively. 



volume of Mo. 



In summary, the Hugoniots of Mo with porosities 
m=1.0, 1.26, 1.83, and 2.31 have been calculated by us- 
ing the first-principles MFP approach. Our results show 
good agreement with the experimental data at m=1.00 
and 1.26. For larger porosities of m=1.83 and 2.31, the 
difference between our results and data points grows with 
the pressure. The Hugoniot temperature of porous Mo 
have also been calculated. For m—1.0, our calculated 
results are remarkably consistent with the previous cal- 
culations. At present, there are no measurement of Hugo- 
niot temperature attainable for porous Mo to confirm our 
theoretical results, and we leave this verification for the 
future shock-wave experiments. 



This work was partially supported by NSFC under 
grant number 10604010 and 10544004. 



1 W. J. Carter, S. P. Marsh, J. N. Fritz, and R. G. McQueen, 
NBS. Spec. Publ. 36, 147 (1971). 

2 H. K. Mao, P. M. Bell, J. W. Shaner, and D. J. Steinberg, 
J. Appl. Phys. 49, 3276 (1978). 

3 J. M. Walsh, M. H. Rice, R. G. Mcqueen, and F. L. Yarger, 
Phys. Rev. 108, 196 (1957). 

4 R. G. McQueen, and S. P. Marsh, J. Appl. Phys. 31, 1253 
(1960). 

5 K. K. Krupnikov, A. A. Bakanova, M. I. Brazhnik, and R. 
F. Trunin, Dokl. Akad. Nauk SSSR 148, 1302 (1963) [in 
Russian] (Sov.Phys. - Dokl. 8, 205 (1963)). 

6 R. G. McQueen, S. P. Marsh, J. W. Taylor, J. N. Fritz, 
and W. J. Carter, The equation of state of solids from 
shock wave studies. - In: High Velocity Impact Phenomena 
edited by R. Kinslow (Academic Press, New York, 1970). 

7 C. E. Ragan, M. G. Silbert, and B. C. Diven, J. Appl. 
Phys. 48, 2860 (1977). 

8 Los Alamos Shock Hugoniot Data, edited by S. P. Marsh 
(University of California, Berkeley, 1980). 

9 L. V. Al'tshuler, A. A. Bakanova, I. P. Dudoladov, E. A. 
Dynin, R. F. Trunin, and B. S. Chekin, Zh. Prikl. Mekh. 
Tekhn. Fiz. 2, 3 (1981) [in Russian] (J. Appl. Mech. Techn. 
Phys. 22, 145 (1981)). 

C. E. Ragan, Phys. Rev. A 25, 3360 (1982). 

1 C. E. Ragan, Phys. Rev. A 29, 1391 (1984). 

2 A. C. Mitchell, W. J. Nellis, J. A. Moriarty, R. A. Heinle, 
N. C. Holmes, R. E. Tipton, and G. W. Repp, J. Appl. 
Phys., 69, 2981 (1991). 

3 R. S. Hixson, and J. N. Fritz, J. Appl. Phys. 71, 1721 
(1992). 

4 R. F. Trunin, N. V. Panov, A. B. Medvedev, Pis'ma Zh. 
Eksp. Teor. Fiz. 62, 572 (1995) [in Russian]. 

5 A. A. Bakanova, I. P. Dudoladov, and Yu. N. Sutulov, Zh. 



Prikl. Mekh. Tekhn. Fiz. 2, 117-122 (1974) [in Russian] (J. 
Appl. Mech. Techn. Phys. 15, 241 (1974)). 
R. F. Trunin, G. V. Simakov, Yu. N. Sutulov, A. B. 
Medvedev, B. D. Rogozkin, and Yu. E. Fedorov, Zh. Eksp. 
Teor. Fiz. 96, 1024 (1989) [in Russian] (Sov. Phys. - JETP 
69, 580-588 (1989)). 

Y. Wang, D. Chen, and X. Zhang, Phys. Rev. Lett. 84, 
3220 (2000). 

Y. Wang, Phys. Rev. B 61, R11863 (2000). 

L. Li and Y. Wang, Phys. Rev. B 63, 245108 (2001). 

Y. Wang and Y. F. Sun, Chin. Phys. Lett. 18, 864 (2001). 

Y. Wang, R. Ahuja, and B. Johansson, J. Phys.: Condens. 

Matter 14, 10895 (2002). 

Y. Wang, R. Ahuja, and B. Johansson, High Press. Res. 
22, 485 (2002). 

Y. Wang, Z.-K. Liu, L.-Q. Chen, L. Burakovsky, D. L. 
Preston, W. Luo, B. Johansson, and R. Ahuja, Phys. Rev. 
B 71,054110 (2005). 

E. Wasserman, L. Stixrude, and R.E. Cohen, Phys. Rev. 
B 53, 8296 (1996). 

J. S. Dugdale and D. K. C. MacDonald, Phys. Rev. 89, 
832 (1953). 

T. Jarlborg, E. G. Moroni, and G. Grimvall, Phys. Rev. B 
55, 1288 (1997). 

P. Blaha, K. Schwarz, G. Madsen, D. Kvasnicka, and J. 
Luitz, WIEN2K, |http://www.wien2k.at71 
P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 
77, 3865 (1996). 

A. C. Mitchell, W. J. Nellis, J. A. Moriarty, R. A. Heinle, 
N. C. Holmes, R. E. Tipton, and G. W. Repp, J. Appl. 
Phys. 69, 2981 (1991). 



